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Abstract 

The collective modes of striped phases in a quantum Hall system are com- 
puted using the time-dependent Hartree-Fock approximation. Uniform stripe 
phases are shown to be unstable to the formation of modulations along the 
stripes, so that within the Hartree-Fock approximation the groundstate is a 
stripe crystal. Such crystalline states are generically gapped at any finite 
wavevector; however, in the quantum Hall system the interactions of mod- 
ulations among different stripes is found to be remarkably weak, leading to 
an infinite collection of collective modes with immeasurably small gaps. The 
resulting long wavelength behavior is derivable from an elastic theory for 
smectic liquid crystals. Collective modes for the phonon branch are com- 
puted throughout the Brillouin zone, as are spin wave and magnetoplasmon 
modes. A soft mode in the phonon spectrum is identified for partial filling 
factors sufficiently far from 1/2, indicating a second order phase transition. 
The modes contain several other signatures that should be experimentally 
observable. 

PACS: 73.40.Hm,73.20.Mf,73.20.Dx 
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I. INTRODUCTION AND SUMMARY OF RESULTS 



Recently, it has been discovered [|I],0 that high quahty two-dimensional electron systems 
in the quantum Hall regime (strong perpendicular magnetic field, low temperature) host 
states with highly anisotropic transport properties. These occur when the filling factor 



V = 2nn£'^ {n is the electron density, i = ^JhcfeB the magnetic length and B the magnetic 
field) is close to half integer with numerator not too small. The strongest effects seem to 
occur for v ~ 9/2, with similar phenomena present at 11/2, 13/2, ... etc. Near these filling 
factors a large asymmetry is observed in the diagonal components of the resistivity tensor 
Pxx and pyy that sets in below approximately 100 mK in GaAs systems. The resistivity ratios 
Pxx/ Pyy may be as large as 3500 0], although the effect is exaggerated by system geometry 
@]. The directions of high/low resistance are clearly correlated with the GaAs crystal axes, 
although the precise mechanism by which they are chosen is at present unknown. For 
states near v = 11/2, 15/2, 19/2, etc, the high/low resistance directions may be rotated by 
application of an in-plane magnetic field [0,0. 

States leading to this anisotropic transport are likely to be related to striped states 
that were found in mean-field studies [^J^ of systems in which several Landau levels are 
filled, and the highest occupied Landau level has a partial filling z/^ in the range 0.35 < 
'^x ^ 0.65. Such ordering has been shown to occur in exact diagonalization studies of 
finite size systems In a seminal theoretical work, the stability of this state to thermal 
and quantum fiuctuations was investigated by Fradkin and Kivelson 0, who pointed out 
a powerful analogy between liquid crystals and quantum Hall stripes. The analogy allows 
a classification of states according to symmetries; these include stripe crystal, smectic, and 
nematic phases. As will be shown below, at zero temperature mean-field theory predicts 
that the stripe crystal is lowest in energy among these. However, it has been argued that 
the smectic state may be stabilized by quantum or thermal [ll| fiuctuations, or both . 



Finite temperature studies |]T2| of a model representing the nematic phase yield impressive 
agreement with experiment of the resistance anisotropy at temperatures that are not too 
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low. Effects of in-plane magnetic fields p3| , p^ have been studied and have provided some 
understanding of the interchange of the high/low resistance directions, although the different 
experimental behavior for v = 13/2, 17/2, ... is still unexplained. 

Beyond transport studies, low- dimensional electron systems may be probed by coupling 
to their collective modes, for example via inelastic light scattering |]TB[ or surface acoustic 
waves |[16|. These collective modes for quantum Hall stripes are the subject of this study. Our 
method will be the time-dependent Hartree-Fock approximation in the form developed by 
Cote and MacDonald fl^. The method requires a static Hartree-Fock groundstate around 
which we can compute excitations. The simplest form for such state is to treat the 
completely filled A^— 1 Landau levels as inert, and form a one dimensional array of alternating 
filled and empty guiding center states in the partially filled Nth Landau level. In this 
approximation, the low energy Hamiltonian for the partially filled level may be mapped to 
the lowest Landau level, with a modified electron-electron interaction. This modification is 
responsible for the low energy of stripe ordering in this system 0. 

We find, however, that uniform stripe states are unstable within the Hartree-Fock ap- 
proximation to formation of modulations along the stripes. The resulting state is essentially 
an ordered array of one-dimensional crystals, i.e., a "stripe crystal" 0. Fig. 1 illustrates the 
charge density for a stripe crystal phase. The amplitude of the density modulation along the 
stripe is small; nevertheless, the energy gained in going from uniform stripes to the stripe 
crystal is considerable. For example, for = 0.5 in the = 3 Landau level, the striped 
phase is found to have energy per particle of —0.279691 in units of e'^/ni (here k is the 
dielectric constant; this will be our unit of energy throughout this paper), while the stripe 
crystal has energy —0.281465. For the parameters of Ref. the energy difference between 
these two states is 112 mK, well above the temperatures for which anisotropic transport is 
observed. Similar results are found at other values of both and u^. 

The energy lowering in forming modulations along the stripes is largely an intrastripe 
effect. For example, one may compute the energy of a stripe crystal with a rectangular unit 
cell rather than the oblique one illustrated in Fig. 1. The modulations of the stripes in this 



state are "in-phase" , requiring an additional Hartree energy. However, due to the weakness 
of the modulations and the long-range nature of the Coulomb interaction, the quantitative 
value of this energy cost is minuscule, of the order lO^^e^/ni ~ 10~^ K. Thus, the chains 
may easily slide past one another. Certainly, at any experimentally attainable temperature, 
the crystal will melt into a series of thermally and quantum disordered one-dimensional 
crystals [|lO| . 



In principle, at zero temperature Hartree-Fock theory predicts the system locks into a 
stripe crystal. The collective modes around this state may be characterized by wavevectors 
k = {k\\,k±), where fcy is the wavevector component parallel to the stripes and k± is the 
perpendicular component. The low energy collective modes are phonons, and in principle are 
gapped everywhere except at k = 0. In practice, because of the small energy scale associated 
with locking, we find nearly gapless modes whenever k\\ =0, independent of k±; the gaps 
are barely resolved by our numerical technique, and are far below currently experimentally 
attainable temperature scales. The low energy collective modes are thus highly reminiscent 
of what is expected for a smectic state 

Fig. 2 illustrates the phonon modes for several values of k± as a function of fcy , computed 
using the time-dependent Hartree-Fock approximation (TDHFA) as described below. Several 
important features are worth noting. 

(1) The modes disperse linearly except for k± = 0, which disperses more slowly. As 
shown below, this is consistent with a harmonic theory of a charged smectic system in a 
magnetic field. The apparent absence of a gap for k± ^ arises because in these collective 
modes the motion of the electrons is parallel to the stripe direction. (This can be seen from 
the eigenvector of the phonon mode from which one can compute the motion of the stripes 
in real space. See Ref. [|I8| for details.) The gap is then controlled by interactions of the 
modulations in different stripes, which is very weak in this system. 

(2) For larger values of A;||, the modes become independent of k±. Physically, this arises 
because the phonon modes are nearly longitudinal for large /cy, involving motion of the 
stripe modulations but no significant motion of the positions of the stripes relative to one 



another |T8|. Since the stripe modulations communicate so weakly, the relative phase of 
motion between stripes has practically no effect on the energy of the mode - hence, no k± 
dependence. 

(3) As might be expected, a gap opens up near k\\ = ±n/b, where b is the distance between 
modulations of a stripe. This leads to a local maximum in the phonon dispersion. In light 
of (2) above, and as may be seen explicitly in Fig. 2, the maximum is extremely fiat along 
the k± direction. As a result, there is a large phonon density of states (DOS) at this energy, 
as illustrated in Fig. 3. Other minima and maxima appear in the phonon dispersion which 
also contribute to structure in the DOS, most notably a double peak at approximately half 
the energy of the n/b peak. Such structures may be observable in inelastic light scattering 
|15|, and their detection would yield optical evidence of stripe ordering in this system. 

(4) A very low energy mode appears along the Brillouin zone boundary at k± = ivr/a 
(a is the separation between the stripes) near k\\ = ±7r/26. As |z/ — 0.5| increases, this mode 
becomes soft (vanishing in energy) just above |z/ — 0.5| = 0.1. This indicates a second order 
phase transition and increased structure in the stripe state as one moves sufficiently away 
from half-filling W^. This may indicate a second order phase transition into a "bubble 



phase" [0 or some precursor of this phase. Alternatively, it may represent a buckling 
instability, in which neighboring maxima within a stripe displace perpendicular to the stripe 
and antiparallel to one another. (Such instabilities are known to occur at Wigner crystal 
edges [^.) The precise motion of the charge in the soft mode is quite complex H^. Work is 
currently underway to determine the precise nature of the groundstate after the instability 
has occured. 

In addition to the low energy phonon modes, the stripe phases support magnetoplasmon 
modes and spin wave modes, and we have explicitly computed them in TDHFA. Fig. 4 illus- 
trates an example of the magnetoplasmon modes appearing as poles of the density response 
function in the first Brillouin zone. The several apparent branches may be understood when 
the structure is compared to analogous modes for a liquid state (no stripes) of the same 
partially filled Landau level, illustrated by the solid lines in the same figure. One may see 



that folding higher order Brillouin zones into the first roughly generates the modes captured 
by the TDHFA. One may thus treat the effect of stripe ordering on these high energy modes 
to a first approximation as that of a periodic potential on an electron gas. A similar effect 
occurs for the spin-wave modes. This does not occur for the phonon modes because no such 
modes exist in the liquid state. 

The presence of several branches of modes near small values of k in principle may be 
detected by optical or surface acoustic wave methods. Such an observation would constitute 
a relatively direct demonstration of striped ordering since it indicates zone-folding effects 
associated with a unidirectional periodic modulation. 

The remainder of this article is organized as follows. In Section |I|, the Hartree-Fock 
method used to generate mean-field states is briefly discussed, and some more details of the 
results are provided. Section |ITT| briefly outlines the method used to obtain collective modes, 
and presents the remainder of our results for collective modes, both in the stripe state and, 
for comparison, in the liquid state. We conclude with a summary (Section [TVI ). There are 
three Appendices. Appendix ^ provides some details of the proper formulation for TDHFA 
in high Landau levels in general and striped states in particular. Analytic expressions for 
collective modes of liquid states for partially filled Landau levels are presented in Appendix 
0. Appendix describes a simple elastic theory demonstrating that the results of the 
TDHFA can be described at long wavelengths by a system with smectic order. 

II. HARTREE-FOCK APPROXIMATION 

In this Section we briefly review the Hartree-Fock approximation (HFA) as developed in 
Ref. [0; some further details are presented in Appendix ^ The fundamental quantities in 
this approach are the operators 

where n, m denote Landau level indices, is the Landau level degeneracy, X are guiding 
center coordinate quantum numbers and a,/3 = ± are spin indices. In the HFA, these 
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quantities are evaluated for a single Slater determinant state, which is accomplished by 
solving the HFA equation of motion for the Green's function [|17 

Gl^^^ (G,r) ^-^Y. {T^n,..Ar)cl^p,x-GyP (0)) exp [-iG^X + ^G.G/V2] , (2) 

X 



with {G} the ensemble of reciprocal lattice vectors of some assumed crystal structure. 
The HFA to the groundstate expectation values of p°'^ (q) are non-zero only for q on the 



reciprocal lattice, and are readily obtained from [1/ 



{pli{G)) = Gt%{G,T = ^-). (3) 
Hartree-Fock energies, electron densities and response functions may be computed from 

{(Pn;^(G)>}. 

For filling factors v = 2N + u, a further simplification/approximation is to project the 
Hamiltonian into the single N — th Landau level, which is formally appropriate when the 
electron-electron interaction scale e'^/ni is much smaller than the cyclotron energy Uc (we 
take h = 1 throughout this paper). While in experimental situations these energy scales 
are comparable, calculations retaining several Landau levels show that, for magnetic 
fields and electron densities relevant to Ref. Landau level mixing lowers the Hartree- 
Fock energy by ~ lO^^e^/zt^ for the striped state. This is sufficiently small to be neglected 
for our present study, and we effectively retain only a single Landau level in our static 



HFA calculations. Assuming also that there is no spin texture in the ground-state, and 



denoting by the index p the partially filled Landau level, we have for < 1 



{Pm (G)) = Sg,o, if m<p] 

(p+(G))^0,(p+(O))=i>, if m=p; 

(Pm(G)) = 0, if m =p; 

{p^{G)) = 0, ifm>p. 



(4) 
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while for v > 1, 



(5) 



{p^ (G)) = 6g,o, if m < p; 

{Pm (G)) = Sg,o, if m = p; 

(p-(G))^0, (p-(0))=l? if m = p; 
^ (p^^,(G)) = 0, if m>p. 

We have defined (G) = p"|° (G) to simphfy the notation. 

With our approximations, the filled levels are inert and cause only a shift of the ground- 
state energy. Up to an unimportant constant, the interaction energy per particle of the 
Hartree-Fock state is then 



^ [Hp, (G) (1 - 5g,o) - X,, (G)] \{p; {G))\' , 



(6) 



G 



where a = + (— ) and z/^. = (z/ — 1) if i/ < 1 {u > 1), and 
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with Jo ix) the Bessel function of order zero and V (q) = 27re^/ g the Fourier transform of the 
electron-electron interaction, for which we use the unscreened Coulomb form. The functions 
[x) are generalized Laguerre polynomials. 

To solve the Hartree-Fock equations, some guess is necessary for the crystal structure 
of the groundstate to specify the set {G}. The simplest structure for the stripes is a one- 
dimensional array with lattice constant a. Writing Cp^o.,x = cx, for a = -|- such states are 
characterized by order parameters 

(8) 



n=— oo 

The density profile of the crystal phase is obtained from the relation 

(^W) = iE(/^p(G)>F,,,(G)e- 



(9) 



G 



where Fp^p{G) is a form factor for electrons in level p (see Appendix A). One can also 
compute a "density" profile corresponding to the guiding centers instead of the real density 
by using 

(^(r))GC = E(/^p(G)>e-^''-'"- (10) 

G 

Such states have been studied for a number of purposes [^^|ll],|l3|,|l4| and provide a 



good first approximation to the Hartree-Fock groundstate at the filling factors of interest. 
However, within the HFA, this state is not stable and cannot be used as a starting point 
for collective mode calculations: the resulting response functions are unphysical. That 
this uniform stripe state is not a minimum of the energy within the space of single Slater 
determinants may be understood as follows. The interaction energy (Eq. (^) for uniform 
stripes may be written as E^f = | Yl'x where ex are the eigenvalues of the Hartree-Fock 
Hamiltonian, and the prime indicates a sum over the Np lowest states, Np being the number 
of particles in the partially filled level. The single-particle spectrum Sx has a well-defined 
Fermi energy Ep with eigenvalues arbitrarily close to it. By introducing a one- dimensional 
modulation along the stripes, a gap is opened at the Fermi energy, the eigenvalues Sx below 
Ep are pushed down, and the total energy is lowered. The resulting state is an array of 
one-dimensional crystals; i.e., a stripe crystal. The collective modes presented below are all 
for such stripe crystal states. 

We conclude this section with some remarks about the results of the HFA. An example 
of the density modulation {n (r)) in a stripe crystal state is presented in Fig. 1. Results for 
other Landau level indices n > 2 and partial fillings Vx ~ 0.5 are qualitatively similar to this. 
Two points are worth mentioning. (1) The amplitude of the density modulations in real space 
are relatively weak, across the stripes and even more so along them. Nevertheless, we will see 
in the collective mode spectra clear signatures of both periodicities. It is interesting to note 
that the weakness of the intrastripe modulations is due mostly to the form factors of the Nth 
Landau level; if one views the "guiding center density" as defined in Eq. (|I^) the intrastripe 
modulations are quite pronounced (cf. flSl). (2) The stripe crystal states studied here break 



particle-hole symmetry; there are separate electron and hole stripe crystal solutions to the 
HFA which at = 1/2 are degenerate. For < (>) 1/2, the electron (hole) crystal is 
lower in energy. 



III. COLLECTIVE MODES IN THE TDHFA 

To obtain the dispersion relation of the collective excitations we compute the matrix of 
response functions 

X^fnlLn, (k + G, k + G', r) = -N^ (Tp";/^^ (k + G, r) p^i^^^ (-k - G', 0)> , 

(11) 

where p (q, t) = p (q, r) — (p (q)) and k is a vector in the first Brillouin zone of the lattice. 
The collective excitations appear as poles of the dynamical response functions and their 
dispersion relation is obtained by tracking these poles for several values of k in the first 
Brillouin zone. Since the order parameters of Eq. (^) were obtained in the Hartree-Fock 
approximation, a conserving approximation for the response functions is obtained in the 
time-dependent Hartree-Fock approximation (TDHFA). In Ref. |T^, it was shown that the 



equation of motion of this matrix of response functions, in the TDHFA, can be written 
schematically as [/ {u + i6) — A] x = B where A and B are matrices that depend on matrix 
elements of the direct and exchange interactions and on the order parameters {{pmi^))} 
only. All response functions can then be obtained by solving numerically an eigenvalue 
equation. In the simplest case (for the intra-Landau level excitation, for example), x consists 
of only one response function and accurate results are easy to obtain. In other cases such 
as for the magnetoplasmon excitations, response functions involving transitions to different 
Landau levels are coupled and the matrix x becomes rapidly very large. Our method is 
thus limited by the size of the matrices x that we can handled numerically. Details of the 
calculation are given in Appendix here we present only the results. For concreteness, 
we focus on a partially filled Landau level of index = 3 and spin a = +, with z7 = 0.45. 
Results for other partial fillings. Landau level indices, and spins are qualitatively similar. 
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To limit the size of the matrix x "we study the collective excitations with uj{k = 0) around 
nuJc ± mg* with m, n = or 1. We assume that uj^ is sufficiently large that coupling 
among excitations near nujc and n'ujc may be ignored iin ^ n' . For comparison, we compute 
the same dispersion relations (when they exist) in the liquid phase i.e. in a homogeneous 
phase with the same filling factor. The dispersions in that case are simply obtained by 
replacing Eq.@ with (^p^ (G = 0)) = u and setting all other order parameters to zero p3[ . 



Many, but not all, of the results we find may be understood in terms liquid-like collective 
modes, whose features have been folded into the first Brillouin zone by the periodicity of 
the striped state. 

There are five types of modes that we consider: (a) n = 0: The phonon mode (present 
in the stripe crystal phase only) appears as a pole of Xnn = Xp^p,p%'^ while the spin-wave 
mode uJsw (k) is a pole of Xo-_ = X^p7p,p'^5 which, according to Larmor's theorem, should 
have 0^51^(0) = g*fiBB. (b) n = 1: There are three magnetoplasmon modes in Xnn that also 
appear in Xa^ = (x^^^^ ~ X^^ ~ X + X ) /4 but with different weight. These 
three magnetoplasmon modes originate from the fact that there are three possible transitions 
with pole around Uc, i.e. (2, -|-) (3, -|-) , (2, — ) — » (3, — ) and (3, -|-) (4, -|-). The Coulomb 
interaction mixes these three modes, with the resulting dispersion branches being quite 
complex even in the liquid phase. A spin-flip mode with 6Sz = +1 appears as a pole of 
X(T+ = X^^^ • The only possible transition is (2, — ) (3, +) and there is correspondingly 
only a single branch in the dispersion. We will refer to this mode as the ujsf+ niode. Finally, 
a pair of spin-flip modes with 6Sz = —I appear as poles of Xa- = ^- These descend 
from transitions of the form (2, -|-) (3, — ) , (3, -|-) — > (4, — ). We will refer to these two 
modes as the usf- modes. 

A. Dispersion relation in the liquid phase 

Figs. 5 and 6 show the dispersion relation of the flve modes for fllling factor u = 0.45 in 
the A^ = 3 Landau level, in the liquid phase. The complex dispersion relations are due in part 
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to the generalized Laguerre polynomial entering in the matrix elements of the Hartree-Fock 
interaction which is responsible for the three minima appearing in all these curves. In these 
figures (and all others that follow), we have substracted the constant term ncUc ± Tng*iiBB. 
Note that two of the magnetoplasmon modes disperses from ujc (as expected from Kohn's 
theorem) while the third one is gapped. The higher energy mode that disperses very rapidly 
is stronger in the density response function Xnn while the lowest-energy mode is stronger in 
Xaz- The middle mode becomes very weak in both response functions as /c — > 0. 

^Prom Fig. 6, we see that the spin wave mode disperses from g*iiBB as expected from 
Larmor's theorem. The inter-Landau-level excitations ujsf+ and ljsf-, however, have their 
gaps cusF+i^) — ~ 9*IJ'bB and cusf- (0) — cUc + g*IJ'BB strongly renormalized by the 
self-energy and vertex corrections. 

B. Phonons 

For the stripe phase, we consider the configuration of Fig. 1 where the electrons on one 
stripe are displaced with respect to the electrons on the other stripes. This stripe crystal 
can be described by an oblique unit cell with one electron or alternatively by a rectangular 
unit cell with two electrons. In the inset of Fig. 2, we show the Brillouin zone of the obhque 
unit cell that extend to k±i = ± 0.44 and to k^i = ± 1.64. 

Unlike the homogeneous liquid phase, the stripe crystal phase can sustain a phonon 
mode. The dispersion relation of this mode is presented in Fig. 2. As discussed in the 
introduction, the most striking feature of the result is the line of nearly gapless modes along 
k± — 0. Generically, for a crystal one expects phonon modes to be gapped everywhere except 
at k = 0. A careful examination of the small limit is consistent with this, although a 
precise determination of the gap is difficult because the mode weights become very small in 
this limit. We estimate the gaps along the k\\ line to be in the range 10~^ — 10~^e^/«;£, which 
is far smaller than any experimentally accessible temperature. Physically, this indicates the 
stripes are free to slide past one another due to thermal fluctuations. 
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In the inset of Fig. 2, we show the dispersion relation along k\\ for several values of 
k±. One sees that the phonons disperse linearly except for = where they disperse 
more slowly. In Appendix C, we show that this is consistent with a harmonic theory of 
a charged smectic system in a magnetic field. Another point worth mentioning is that 
for larger values of k\\ the dispersion in k± becomes almost independent of k^. By direct 
examination of the charge motion in several such collective modes, we have found that 
this arises because the phonon modes are nearly longitudinal for large k\\; they do not 
involve significant motion of the positions of the stripes relative to one another. Since the 
stripe modulations communicate so weakly, the relative phase of motion between stripes has 
practically no effect on the energy of the mode. As discussed in the Introduction, this results 
in resonances in the collective mode density of states that might be observed in inelastic 
light scattering. A particularly strong such resonance occurs due to the additional flatness 
of the dispersion near the Brillouin zone boundary along the direction of the stripes (see 
Fig. 2) where a gap opens up separating the "acoustic" from the "optical" modes. 

Finally, as discussed in the Introduction, a soft mode appears that indicates an instability 
of the modulated stripe state studied here for z/ just below 0.40, suggesting at these lower 
fillings that the correct groundstate will have more structure. We note that this instability 
indicates a second-order transition, in contrast to the first-order transition found near u ~ 
0.36 between stripe and "bubble" states studied in Ref. [Q. The result may indicate that 
a precursor of the bubble phase develops within the stripe phase, perhaps in which the 
bubbles are elongated rather than circular. It is also possible that the stripes have a buckling 



instability in analogy with similar behavior previously noted for Wigner crystal edges |20 



C. Higher Energy Modes 

Unlike the phonon mode, the four other excitations that we consider also exist in the 
liquid phase. To understand the effect of the stripes, we plot the dispersion relations ob- 
tained in the stripes and liquid phases together. A few comments on these results are 
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in order before we present them. For the hquid, we fold the modes in the first Brillouin 
zone of the stripe crystal and keep the lowest-energy branches {u; (k + G)}. Along the di- 
rection perpendicular to the stripes, the lowest-energy branches correspond mostly to the 
functions uj {k±_ -\- nG±, k\\ — O) with n — 0, ±1, ±2, ... In the direction of the stripes, they 
correspond mostly to the curves cu (k± — nG±, k\\j . In this case, however, the curves with 
n = ±1, ±2, ±3, ... have the same energies in the liquid; i.e., are degenerate. (The thick lines 
in the figures represent n — which is not degenerate) . This degeneracy is sometimes lifted 
in the stripe phase. Note that this Brillouin zone folding of the liquid dispersions introduces 
a large number of branches. It is not possible to track all the corresponding poles in the 
stripe phase. We thus sometimes show only a small subset of these modes corresponding 
to low-energy excitations. Because we keep only the most intense poles and because the 
relative intensity changes as k spans the Brillouin zone, the dispersions sometimes appear 
discontinuous for the stripe phase; this is because the mode weights fall below our threshold 
for plotting them. Note that the zone-folding effects lead to the presence of several branches 
near small values of k. In principle, these excitations could be detected by optical or sur- 
face acoustic wave methods and would thus represent a direct demonstration of the stripe 
ordering. 



1. Spin Waves 

For GaAs systems, under most circumstances g*iiBB « hwc, so that spin waves are 
the lowest energy modes after the phonons. Figs. 7 and 8 show the dispersion relation 
obtained for the spin waves. The most striking difference between the zone-folded liquid 
results and the spin waves of the stripe state is a dramatic anisotropy in the gap opening 
at the Brillouin zone boundary. This gap is much larger at the boundary for large k± 
than the corresponding one for large k\\. Certainly, part of the explanation is that the 
density modulations responsible for the latter is much smaller than that of the former. 
However, the gap at large k± is much larger than, for example, the corresponding gap in 
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the magnetoplasmons, discussed below. This strong many-body effect may be related to 
electrons at the stripe edges having a small local spin stiffness relative to those in the liquid 
state or in the center of the stripes. In any case, this many-body effect results in a branch 
of spin waves with surprisingly low energy. 

2. Magnetoplasmons 

Figs. 4 and 9 show the dispersion relation obtained in the stripe phase for the magneto- 
plasmon modes. To capture all the three branches, we show poles obtained from both the 
density and spin response functions Xnn and ■ Because the three corresponding branches 
in the liquid are almost flat at large wave vector, the folding of the modes in the first Bril- 
louin zone introduces many branches at small energy. It is quite clear, however, that the 
dispersion obtained in the stripe phase follows closely that of the liquid, with small gaps at 
the Brillouin zone edges and some lifting of degeneracy in the k\\ direction. 

3. Spin Flip Excitations 

Some of the spin flip excitations seem to follow behavior reminiscent of our results for the 
magnetoplasmons, closely following the liquid results, whereas others undergo strong many- 
body renormalizations, as we found for spin waves. Fig. 10 shows the dispersion for the cusf- 
mode. (In these figures, we show only the low-energy excitations because the liquid modes 
become very complicated at higher energies.) One can see the direct correpondance between 
the liquid phase dispersions and those of the stripe states. As for the magnetoplasmons, 
these differ by gap openings and the lifting of degeneracies. For the u)sf+ rnode, as for 
the spin wave mode, the dispersion relation of the lowest-energy branches are quite close to 
the corresponding liquid result. For higher branches, however, the stripes ordering lead to 
important changes as can be seen in Fig. 11. The difference between the behaviors of ojsf+ 
and u!sF- is very likely related to the presence of two branches of the former in the liquid 
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state which may be mixed by the stripe ordering, whereas only a single branch exists in the 
latter. 



IV. SUMMARY 

In this work we have studied collective mode of stripe states for quantum Hall systems. 
The lowest energy states are phonons, with a line in the Brillouin zone of extremely low 
energy states, making the resulting low-energy physics of this system that of a charged, two- 
dimensional smectic in a magnetic field. We also found signatures in the phonon density 
of states indicative of stripe ordering that should be detectable in inelastic light scattering, 
and a soft mode that indicates an instability of the stripe state for partial fillings sufficiently 
far from 1/2. Results for spin waves, magnetoplasmon, and spin-flip excitations were also 
presented, which in a first approximation could be understood in terms of zone-folding of 
corresponding excitations for the liquid state. Some of these, however, underwent strong 
renormalizations due to electron-electron interactions; in particular, we found a surprisingly 
low-energy branch in the spin wave spectrum due to this effect. 

The form of the low energy physics of this system has important consequences for quan- 
tum fluctuation effects on the stripe crystal state, particularly the stability of the crystal as 
well as pinning by disorder. Some of this has been discussed previously |10|; a more detailed 
study will be presented in future work. 
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APPENDIX A: DETAILS OF THE TDHFA 



In this Appendix we discuss the proper formulation of the TDHFA in high Landau levels. 



The basic approach follows that of Ref. |17]; however, there are important details involved 



m 



in computing inter-Landau level excitations that have been treated incorrectly |T7.24 
the literature leading to results that do not correctly include the exchange self-energy 
corrections to these excitations. We present here a correct formulation of the TDHFA that 
avoids such errors and respects Kohn's theorem 



1. Static Hartree-Fock Approximation 

We begin by briefly reviewing the relevant equations for HFA that will be needed in our 



formulation of the TDHFA; details may be found in Ref. |17]. Our model HF Hamiltonian 
is 

Hhf = N^Y1 ^"--Z^"-- + E E (G) Pn (G) , (Al) 

n,a n,a G 

where 

en,a = {n + 1/2) uJc - ag*fiBB/2. (A2) 
The Hartree-Fock effective potential (G) is given by 

U: (G) = E E ("^' ^' ^' ^5 - 5^,pX (m, n, n, m; G)] {pi (-G)> (A3) 

m f3 

= E E [^-'- (G) - 5«,^A„,„ (G)] {pi (-G)> . 

m 13 

For completeness, we give here the form of the Hartree and Fock interactions that enter into 
the calculation of the self-energy corrections to the collective excitations: 

br n < m). 



nil qi V 2 



X™,n(q) = (^) v^(£t) dxx'^^--^e-^' [L^^ {x')]' Jo (v^xg^) (f. 
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For n > m, we use Xn,ni (q) = -^m,n (q)- 

The effective interactions appearing in Eg. (|A3[) are a subset of tlie more general form 



H (m, n2, na, n^] q) = ( (^77-^ ) ^ (q) ^ni.na (q) Fn^^n^ (-q) , (A5) 



X(ni,n2,n3,n4;q) = j (q') (q') F.3,n. (-q') e-^''^^''^ (A6) 

that we need to derive the TDHFA. We take V (q) = 2Txe^/q. (We use the two-dimensional 
cross product as a short form for q x G =qxGy — QyGx)- These interactions contain the form 
factors 



1/2 / / ■ n\\ n—m 



m\\ ' ( {-Qy + iqj 



-.2 02 



i;r"' ( ^ I ■ (A7) 



for m < n, where (x) is the generalized Laguerre polynomial. Note that Fn^m (q) = 
[Fm,n (~q)]* • We remark that the effective potential in any Landau level n,a depends on 
the occupation of the other levels, as does the energy of the electrons in that level. This 
self-energy shift differs from one level to another, and makes an important contribution to 
the energy of inter-Landau level excitations. 

The single particle Green's function of Eq.(|^) obeys, under the Hamiltonian of Eq. {"Kl), 
the equation of motion 

[z^„ - (e^ - /i)] (G,^^„) - ^ W: (G - G') iG',tiUn) = 5g,o, (A8) 

G' 

where /x is the chemical potential and 

(G - G') = (G' - G) e^Gx<^'^V2_ (^9) 
Eq. {\KBi) can be solved numerically to compute the densities (p^ (G)) as explained in Ref. 
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2. Time-Dependent Hartree-Fock Approximation 



The two-particle Green's functions are defined by Eq. (0). In the TDHFA, they obey 
an equation of motion that we write as fl] 



+ (e^,,^ - xi?,nf;ll (k + G, k + G', (AlO) 

+ E bG,G" (k) (G" - G) - 7G,G" (k) f/4 (G" - G)] (k + G", k + G', fi„) 

G" 

= K.,nA.,nA,sSp,, [7g,g' (k) (p°, (G - G')> - 7G,G' (k) {pi (G - G')>] , 
x;^&l3,n, (k + G, k + G', Qr.) = xKSl, (k + G, k + G', (All) 

- E EE^Sf;rne(k + G,k + G'',a„) 

n5,---ng r?i/ G" 

X X {nj, ne, ns, ng; k + G") X^fna.n, (k + G", k + G', , 

and 

XnSl.n. (k + G, k + G', Qn) = Xnfnl'Ln, (k + G, k + G', Q.^) (Al2) 

+ E EE^™,n.(k+G,k + G'',fi„) 

n5,---ng »7!/ G" 

X H (n5, ne, nr, n^; k + G") x'^^::f,n,,n, (k + G", k + G', , 
where f2„ is a Boson Matsubara frequency and 

7G,G'(k)=e*('^+^)x('^+^')^'/2_ ^A13) 

Eqs.(AlO - A12) are equivalent to the result of summing ladder and bubble diagrams in a 
perturbative expansion of x- Note that the only information required in these equations 
is the groundstate density < p°(G) >. The equations couple together an infinite set of 



response functions; as discussed in Ref . ||T^ , when truncated appropriately they may be cast 
in a matrix form for numerical solution. In the next few sections, we describe truncations 
and simplifications that are appropriate for computing various collective modes. 
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3. Equation of motion for x^^"^ 

It follows from Eq.(AlO) that the only non zero x^^"* niust be of the form Xn}n,m,n'" ■ Writ- 
ten in matrix notation (with the reciprocal lattice vectors G, G' being the matrix indices), 
the equation of motion for x^^^ is then 

[znj - F:;i (k)] (k, u) = BZ:i (k) , (A14) 

where 

[f:;^ (k)] ^ ^ {e^,p - 5g,g' - K (G' " G) 7g,g' (k) + f/^ (G' - G) 7g,g' (k) 

(A15) 

and 

[Bui (k)] ^ = 7g,g' (k) {Pl (G - G')) - 7g,g' (k) {pi (G - G')> . 

(A16) 

The size of these matrices depends on the number of reciprocal lattice vectors that are kept 
in the numerical calculation. 

4. Equation of motion for x 

Since x''"'* is of the form Xn)rn,m,n", Eq.(All) for x can be simplified to 

X-nfnl'Ln, (k, ^) = xS^^Kn", (k, ^) 5„,„,5„3,n.5«,55A, (A17) 

_ (0),a,/3,/3,a /i^ \ ST^ I v- ~a,/3,7,<5 I'l^ 

Ani,n2,n2,ni K'^^ ^ ) / ^ [^»-n5,ni,n2,Ti6 \'^) A.n5,n6,n3,n4 K'^^^Jj ' 

where the matrix 

[^.„„2,.3,.4 (k)] G,G' = ^ (^1' ns, n,- k + G) 5g,g'- (A18) 

lYiom. Eq.( |A17| ), it is clear that, in the spin indices, x must be of the form x"'^'^'" and 
that Xnf,n2^n3,n4 7^ 0' ^uly if Xnln^^'^nf^'ni 7^ 0. Siucc wc are working in the strong magnetic field 
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limit {uJc » e^/ni), we will assume that a response function with poles around nUc is only 
coupled to other response functions poles near the same frequency |2^. Thus, we truncate 
our equations by including coupling among response functions of the form Xni+m,n2+m,n3,n4,- 
for different values of m. We remark here that coupling all response functions with pole 
around the same frequency is essential to recover Kohn's theorem for the cyclotron mode. 
Eq. (|A17|) is now simplified to 



Xni+m,n2+m,n2 



my^J^) Xni,n2,n2,ni y^^^ ) "rn,0 Xni 



(0),a,/3,/3,a /i x 

+m,n2+m,n2+m,ni+m ^ ) 



(A19) 



X 



^ ] ^^ni+m' ,ni+m,n2+m,n2+m' (k) Xni+m' ,n2+m' ,n2,ni (k, 



Because we consider only the special case where all Landau levels below p are completely 
filled and p is partially filled, Eq. ( [Al(j| ) implies that 

x!p}-m,p+m,p+m,p+m 7^ ^ Ouly if 777, = 0, 

(0),a,/3,/3,a , „ , ;f _ _i n 

Xp+m,p+m+l,p+m+l,p+m T~ ^^^^J '"' — ^I'^i 

x\,+m,p+m+2,p+m+2,p+m 7^ Ouly if TTT, = —2, —1, 0, 



(A20) 



and so one. x will thus be coupled to one, two, three or more other x's depending on the 
value of m and also on the number of levels filled below p. For example, we only need to 
consider the response function for the intra-Landau-level excitation. Its equation of 

motion is thus simply 



Xpitp (k, ^) = (k, ^) - (k, ^) Xf^ (k) x^t^'p" (k, ^) , 



where we have defined the diagonal matrix 



(A21) 



(k)] ^^,=X (n, n, n, n; k + G) 5g,g'. 



(A22) 



With Eq.([mD, Eq.([MTD becomes 



[^fi„/ - F^f (k) + (k) X(°) (k)] Xptf (k, u) = B^f (k) 



(A23) 



21 



For m — 1 (inter-Landau-level excitations), there are only four non zero x with poles 
around +cUc- To deal with this case, it is helpful to define the block matrices (which we 
denote by the symbol x to distinguish it from the simple matrix x) 



X 



and 



-^p,p+l,p+l,p -^p,p+l,p,p— 1 

~a,l3,l3,a ~a,l3,l3,a 
A.p-l,p,p+l,p Ap-l,p,p,p-l 



(A24) 



1 = 



I 
I 



(A25) 



In terms of these matrices, Eq.(All) simplifies to 



where 



and 



with 



iO„/-F"'^(k) + B"'^ {k)Xp{k) 



X (k, u)^ B (k) , 



(k) 



i^p-llp (k) 



B""'^ fk) = 



Bpt^. (k) 



B;'X^ (k) 



(A26) 



(A27) 



(A28) 



(A29) 



W(k)]G,G' 

[^f(k)]c,c' 

[^^^ (k)] 



G,G' 



X (n, n, n + 1, n + 1; k + G) 5g,g'', 
X {n,n — l,n,n + G) Sg,g'] 

X {n — l,n — l,n,n;k + G) 5g,g'', 
X (n — 1, n, n + 1, n; k + G) Sg,g'- 



(A30) 
(A31) 
(A32) 
(ASS) 
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The solutions to Eg. ( |A26|) can be used to compute both density and spin-flip response 
functions (e.g., Xp.p+w,p)- 

In principle, we can deal in the same manner with excitations around 2ujc. These would 
involve solving a 3 x 3 block matrix in x (depending upon the number of filled levels below 
p). Since each block in these matrices is itself a matrix whose size depends on the number of 
reciprocal lattice vectors that we keep in the calculation, solving for higher-energy excitations 
becomes difficult numerically. We will thus be satisfied here with the solution for intra and 
inter-Landau-level response functions with poles around zero or Uc (shifted, of course, by 
the Zeeman energy if spin fiip excitations are considered). This includes the important case 
of phonons and spin-wave excitations in the partially filled level, as well as the cyclotron 
modes around Uc and spin-fiip modes from or to the partially filled Landau level. 



5. Equation of motion for x 

The full response function x is computed by including the Hartree vertex corrections 
which, from Eq.(Al2), gives 

XnSla.n, (k, ^) = X:f£%,,n, ^) (A34) 
+ ^Q,/3 Xnun^'^5,n6 (^j ^) ,n6,n7,n8 (k) Xnr^nf,n3,n4 (k; ^) ! 

where 

[H„^,r.„n„n, (k)] ^ = H {ui , n2 , , u^; k + G) 6g ,G' ■ (A35) 

To simplify this equation, we will again consider separately the case of intra and inter- 
Landau-level excitations. 

For intra-Landau-level excitations, we write 



/\,n,n,m,m'- 



E\ ^ ~a,a,a,a tt 
/ J Xn,n,n5,n(i-'^n5,n(i,nj,ns 



(k) x"'"'''' 



nT,ns,m,m 



{k,u;) 



(A36) 
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Using our approximation of no coupling between excitations of different noOc, one may show 
that 



[^^ni - {f;:; (k) + s;^- (k) [h^^ (k) - (k)] )] xZZ (k, ^) = ^.T (k) • 



(A37) 



where a is the spin index of the partially filled level and 

(k)] ^ = H (n, n, n,n-k+ G) 6g,g'- (A38) 

The poles of Xp'p'p'p contain the phonon mode. For spin waves, there is no Hartree vertex 
corrections and the relevant response function obeys 

[tQJ - (k) - B;:- (k) (k))] 4^;^+ (k, = B^f (k) . (A39) 

The most complex situation is that of inter-Landau-level density modes (cyclotron 
modes). In this case, we need to consider the coupling between Landau levels as well as 
between spins. From Eq. (A34) and Eq. ( [A2(j| ), we get 

[tnj - (k) + 5"'" (k) Vp (k)] T'"'^'" (k, uj) 



(A40) 



i?"'"(k)5„,;3 + i?"'° (k) (k) 



where 



and 



a,a,f3,f3 a,a,f3,f3 

^p,p+l,p+l,p Xp,p+l,p,p— 1 

a,a,f3,f3 a,a,l3,f3 

Xp-i,p,p+i,p Xp-i,p,p,p-i 



^(1) ^(4) 



(A41) 



7t{2) o-(3) 



(A42) 



with 



G,G' 



[^f (k)]o,c' 
W^(k)]o,c' 



G.G' 



H {n+ 1, n, n, n + 1; k + G) (5g,g'; 
i7 (n, n - 1, n, n + 1; k + G) 5g,g'; 
i7 (n, n — 1, n — 1, ri; k + G) 5g,g'; 
(n + 1, n, 77, — 1, n; k + G) 5g,g'- 



(A43) 
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Eq. (A40) can be written in a more transparent form by defining the block matrices 



X = 



,+ 



Xp-l,p,p+l,p Xp-l,p,p,p-l Ap-l,p,p+l,p Ap-l,p,p,p-l 
Ap,p+i,p+i,p Ap,p+i,p,p— 1 Ap,p+i,p+i,p /Cp,p+l,p,p— 1 
Xp-i,p,p+i,p Xp-i,p,p,p-i Xp-i,p,p+i,p Xp-i,p,p,p-i 



(A44) 



B 



B ' 
B~' 



B. 



p>p+i 
B, 








(1) 



Hp — Xp 



H 



(2) 



1^(1) o-(4) 

vVp ilp 

1^(2) rr(3) 

y\p Hp 



(4) 



H 



(1) 



(2) 



(4) 



(3) 



p-l,p 










B, 







'p-l,p 



H 



(1) 



if, 



(4) 



H 



(2) 



(3) 



-Hp ^ ~ Xp 



(1) H-(4) 



Hp — Xp 



H 



(3) 



(4) 



o-(2) _ ^(2) „(3) _ ^(3) 



(A45) 



(A46) 



(A47) 



/ 

7 



J 
0/00 

7 

/ 



(A48) 







F' 





i^: tp (k) 

OF, 










p,p+i (k) 

i7-r,p (k) 



(A49) 
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The equation of motion then takes the form 

'tQj - f (k) - 1 (k) % (k)] f (k,.;) = 1 (k) (A50) 

We solve this matrix equation numerically to obtain response functions whose poles give the 
magnetoplasmon and inter-Landau-level spin-flip excitations. 

APPENDIX B: DISPERSION RELATIONS IN THE LIQUID PHASE 

The equations of the previous section can be drastically simplified in the homogeneous 
or liquid phase since then (p" (G)) = i^^co or {p^ (G)) = 0. For example, for < z/ < 2, the 
dispersion relations for various modes can be computed analytically, and one may show they 
reproduce the results of Kallin and Halperin . Larger filling factors are more complicated 



as they involve several different particle-hole excitations; the coupling among these has 
not been treated correctly in previous studies [^]. As concrete examples of the present 
method, we compute the collective modes for the liquid state at different filling factors 
and, in particular, for u = 6.45 which corresponds to the filling factor of the stripe crystal 
considered in this paper. 

a. Liquid phase with < < 1 

As an application of the above formalism, we consider here the simple case of < z/ < 1. 
If the lowest Landau level is partially occupied with up spins, then (p^}" (G)) = 1^60,0- AH 
matrices are diagonal and so k + G — q which is not restricted to the first Brillouin zone. 

Since BqI^ (q) = if follows that there can be no phonon mode. Moreover, since 
Bq I^ (q) = z/, we have from Eq. ( |A39| ) 



uJsw (q) = 9*f^BB + V Xo,o (0) - xf (q) . (Bl) 



Since xf^ (q = 0) = Xo,o (0) , it follows that ^sw (0) = g*PBB as required by Larmor's 
theorem. 
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For the density mode Bq'^ (q) = u, and the dispersion is 



l^nn (q) = + 



Xo,o (0)-Xi,o(0) + iff) (q)-Xr (q; 



'(0) 



(B2) 



In this equation z/Xq.o (0) is the self-energy lost by the electron leaving level n = while 
z/Xi (0) is the self-energy gained in the new level n = 1. Because h[°^ (0) = and Xo,o (0) - 
Xi,o(0) - Xf ^ (0) = 0, it follows that cu^yc. (0) = ujc as required by Kohn's theorem. We 
remark that these two results are identical to those of Ref. 

For the inter-Landau-level spin flip excitation, Bq^^ (q) = and the dispersion is 

usF (q) = + g*f^BB + V [Xo,o (0) - xf^ (q) 



Notice that 



(B3) 

Xo,o (0) - Xf^ (0) > so that the self-energy orrection introduce 

a positive shift in the dispersion relation contrary to the result in Ref. This was first 
noticed in Ref. [^. Similar problems with the inclusion of the self-energy terms appear 



in the higher-energy modes as well in Ref. |2J]. Apart from this discrepancy, our results 
reproduces correctly the dispersion relation of the higher-energy modes of the liquid phase. 



h. Liquid "phase with 6 < < 7 

This is the case we consider in the stripe phase. It is thus interesting to compare the 
dispersion relations obtained there with the corresponding ones in the liquid phase. We 
assume that the partially filled level is {p = 3, +), so that (pg (G)) = z/5g,o and (p^ (G)) = 
5g,o for Tn < 3. There is again, of course, no phonon mode. 

For spin wave, -63^3" (q) = and the dispersion, from Eq. (|A39 ) is simply 



u^sw (q) = g*l-iBB + V [X3,3 (0) - Xf (q)J , (B4) 

with ujsw (0) = g*fiBB as required. 

For spin-flip excitations with 6Sz = +1, we must look at X2 3*3^'" which, from Eq. (|A26|) 



is coupled to X3^4|3^2' • Solving this system of equation, we rapidly obtain that Xs,'^^!^' = 
so that the dispersion relation is obtained from X2 s^s't " only. This makes sense, since the 
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transition (2, — ) (3, +) is not coupled to any other in the situation we consider. We find 
then 

LUSF+ (k) = ujc- 9*I^bB + J:^;^^ -(!-/>) (k) , (B5) 

where 

^2,3^ = X2,0 (0) + X2,l (0) + X2,2 (0) - X3,o (0) - Xg,! (0) - Xs,2 (0) - uX^^S (0) . 

(B6) 

The exchange and vertex corrections introduce a downward shift in usf+{0) from ojc—g* [j^bB 
(see Fig. 6). 

For spin-flip excitations with 5Sz = —1, the transition (2,+) — > (3,—) is coupled to 
(3,+) — s> (4,—) and we must solve Eg . (|A26|) that couples X3,'4~4:^'~^ to X2^374,3'~'^- There is 
correspondingly two such spin-flip modes, with dispersion given by 

usF- (k) =uj, + g*fXBB+^ {At;,- (k) + Aj- (k)) (B7) 
± y (Ajs- (k) + Ajr (k))' - 4 [Ajr (k) Ajr (k) - l^X^'^ (k) (k)' 
where 

A+3- (k) ^ S+3- - z7X« (k) , (B8) 
At:, (k) = Sj,- - Xf (k) , (B9) 

with 

Sjs" = X2,o (0) + X2,i (0) + X2,2 (0) + Z>X2,3 (0) - Xg.o (0) - Xg,! (0) - X3,2 (0) , (BIO) 
= ^3,0 (0) + Xs,! (0) + X3,2 (0) + uXs,3 (0) - X4,o (0) - X4,i (0) - X4,2 (0) . (Bll) 

In this case, the shift is positive in both modes (see Fig. 6). 

For the density modes, three excitations are coupled: (2, +) — > (3, +) , (2, — ) (3, — ) 
and (3, +) (4, +). Since -B^4~ (k) = 0, X3,474,3'~ (k, i^) = and the 4x4 block matrix in 
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Eg. ( |A44| ) reduces to a 3 x 3 block matrix. The three collective modes are found from the 



determinant of 



[uj + i5)I -F{k)-B (k) Vp (k) in Eq.(^) i.e. from 



- ^c) - At:: (k) (k) - (k)) -uHP (k) 

(1 - u) (i/f (k) - Xf (k) 



where 



and 



(UJ — iUr 



Ks^ (k) 



hP (k) 



K: (k) 

A^Y (k) 



A.r (k) 



^3,4 
^2,3 



^2,3 



i/f ) (k) 



(k) 

{uj - cuc) - Ao r (k) 



= 0, 
(B12) 



i7 (ijf (k) - xi^) (k)^ 

(1 - u) (//f ) (k) - Xf ) (k)) 



i/i^) (k) - Xr^ (k) 



(3) 



(B13) 
(B14) 
(B15) 



^3^4^ = ^3^4 ~ ^^-^4,3 (0) , 
Sjs^ = - Z7X3,3 (0) , 

S2,r = ^2,0 (0) + X2,i (0) + X2,2 (0) - X3,o (0) - Xg,! (0) - X3,2 (0) 

These collectives modes are represented in Fig. 5. 

APPENDIX C: HARMONIC THEORY 



(B16) 
(BIT) 
(B18) 



As discussed in the text, to an excellent approximation the phonon mode frequencies 
computed in the TDHFA for the stripe phase disperse linearly from fcy = 0, with a slope 
that vanishes at = 0. In this Appendix we demonstrate that this behavior is consistent 
with a harmonic theory for a two-dimensional charged smectic system in a magnetic field. 
Defining the direction parallel to the stripe as the x direction, the simplest long-wavelength 
harmonic potential one can write down might be 



K 



dx 



duy 

dy 



(CI) 
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In the above equation, u is a displacement field for the stripes, the first term represents an 
elastic contribution for longitudinal compression of the stripes, and the second arises from 
interstripe repulsion. Collective modes are most easily computed in terms of the Fourier 
transformed displacements, 

. . 1 



rfVe*'i "u (r). 



(C2) 



where A is the area of the system. To compute the collective modes in a single Landau level, 
one may impose the commutation relations |^ [u^ (qi) {^2)] = '*^^'^qi,-q2- The equation 
of motion = \u^, V~\ , where fi = x,y, after Fourier transform with respect to time may 
be written in the form 

/ 



\ 


(q) ^ 
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(C3) 



Dyx (q) Dyy (q) 



For a system with inversion symmetry, D^y (q) = Dyx (q), and the eigenvalues of Eq. (|C3| 
(i.e., the collective mode frequencies) take the general form 



q) = ±tjDxx (q) Dyy (q) - \Dxy (q) P. 



(C4) 



, ^^^^ = K^qj^, D,j.y = 0, so (q) oc QyQx- This has the correct behavior of linear 



ForEq. (fn\),D 

dispersion with respect to qx, with a slope that vanishes as qy —>■ 0. However, the model has 
the incorrect behavior of containing zero modes along both the qx and qy axes p8| 



The key missing ingredient in Eq. ( |(J1D is that the restoring force for motion perpendicular 
to the stripes comes only from interstripe repulsion. We expect, however, that individual 
stripes resist bending, as in a smectic system. Thus, one should add a curvature term to 
the energy, which in this case we take to be 



1 f f (Puy^ ^ 



(C5) 



Adding this to the resulting Dyy (q) is modified to i^yq^ + i^h<it- The collective mode 



frequencies then take the form u; (q) oc qxy<ly+ (^^^ it- This has the correct behavior 
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that collective modes are gapped except for = 0. A further refinement necessary to 
correctly describe the long-wavelength physics of this system is the addition of the Coulomb 
interaction. This may be simply modeled by adding a term of the form 

2 



(C6) 



Vcoul 

■I ^ q 
q 

with Kc = 27re^/Ka^, where in this last expression k is the dielectric constant of the host 
material, and Oc is the area per electron in the groundstate. Finally, symmetry also allows 
the addition of a term of the form 

Taking our potential energy to be ^ + Hend + Vcoui + Vxy, one may easily compute D^^^, 
and use Eg . (|C4|) to show uj (q) oc for qy > 0, and uj (q) oc ql^'^ for qy = 0. The sublinear 



behavior of the collective mode at = is consistent with the results of the TDHFA. 
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Figures Captions 

Fig.l Electron density for a stripe state with v = 6.45. The separation between the stripes 
is a = 7.16£ and the period of the modulations along the stripes is 6 = 1.95^. The mod- 
ulations on two adjacent stripes are displaced by 6/2. This pattern can be described 
by a primitive unit cell with lattice vectors Ri = (a, 6/2) and R2 = (0, b) . 

Fig. 2 Phonon dispersion relation of the stripe phase with 1/ = 6.45. Left inset: phonon 
dispersion along the stripes for different values of the wave vector k±. The local 
minimum has a frequency of a; 0.02 (e'^/ni) and becomes soft as fiUing factor is 
decreased. Right inset: Brillouin zone for the primitive unit cell described in Fig. 1. 

Fig. 3 Phonon density of states. (The small oscillations are numerical artifacts). 

Fig. 4 Dispersion relation a; — of the three branches of the magnetoplasmon mode in the 
stripe phase along the direction perpendicular to the stripes in the density pattern 
of Fig. 1. The corresponding dispersions in the liquid phase (see Fig. 5) have been 
folded in the first Brillouin zone and are represented by full lines with the heavy lines 
indicating parts of these dispersions that lie in the first Brillouin zone. The filled circles 
represent frequencies obtained from Xnn while the empty squares represent frequencies 
obtained from Xaz ■ 

Fig. 5 Dispersion relations cu — cUc oi the three magnetoplasmon modes of the liquid phase 
with v — 6.45. 

Fig. 6 Dispersion relations of the spin wave and spin flip modes for v — 6.45 in the liquid 
phase. The dispersion usw " g*I^BB of the intra-Landau-level spin wave mode is 
represented by the dashed line. The dispersions ujsf- — {^^c + 9*IJ'bB) of the two 
branches of the inter-Landau- level spin flip mode with 5Sz — —1 are represented by 
the full lines. The dispersion ljsf+ — {^c ~ g*l-'-BB) of the inter-Landau-level spin flip 
mode with SSz — +1 is represented by the dot-dashed line. 
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Fig. 7 Dispersion relation usw — g*I^BB of the spin wave mode in the stripe phase along 
the direction perpendicular to the stripes in the density pattern of Fig. 1. The 
corresponding dispersion in the liquid phase (see Fig. 6) has been folded in the first 
Brillouin zone and is represented by full lines with the heavy line indicating parts of 
these dispersions that lie in the first Brillouin zone. 

Fig. 8 Dispersion relation cusw ~ 9*IJ'bB of the spin wave mode in the stripe phase along the 
direction parallel to the stripes in the density pattern of Fig. 1. The corresponding 
dispersion in the liquid phase (see Fig. 6) has been folded in the first Brillouin zone 
and is represented by full lines with the heavy line indicating parts of these dispersions 
that lie in the first Brillouin zone. 

Fig. 9 Dispersion relation a; — a;c of the three branches of the magnetoplasmon mode in the 
stripe phase along the direction parallel to the stripes in the density pattern of Fig. 
1. The corresponding dispersions in the hquid phase (see Fig. 5) have been folded in 
the first Brillouin zone and are represented by full lines with the heavy lines indicating 
parts of these dispersions that lie in the first Brillouin zone. The filled circles represent 
frequencies obtained from Xnn while the empty squares represent frequencies obtained 
from Xa,- 

Fig. 10 Dispersion relations cusf- — (<^c + 9*IJ'bB) of the spin flip mode in the stripe phase 
along the direction perpendicular to the stripes in the density pattern of Fig. 1. The 
corresponding dispersion in the liquid phase (Fig. 6) has been folded in the first 
Brillouin zone and is represented by full lines with the heavy line indicating parts of 
these dispersions that lie in the first Brillouin zone. 

Fig. 11 Dispersion relations u!sf+ — {^c ~ g*I^BB) of the spin flip mode in the stripe phase along 
the direction parallel to the stripes in the density pattern of Fig. 1. The corresponding 
dispersion in the liquid phase (see Fig. 6) has been folded in the first Brillouin zone 
and is represented by full lines with the heavy line indicating parts of these dispersions 
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that lie in the first Brillouin zone. 
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